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ABSTRACT 

The objective of this study is to investigate Vertical 
Steam Condenser design optimization, in an effort to develop 
condensers with minimum weight and volume. Candidate vertical 
condenser tubes considered were of two types: (1) external 

axial flutes with smooth internal tubes, and (2) doubly en- 
hanced tubes that include internal (coolant side) heat transfer 
augmentation. A computer program (VERTCON-2) was developed to 
be used as a preliminary design tool for selection of the 
optimal condenser that meets the specified geometrical, and 
operating parameter limits. 

The computer model considers condensation in the 
presence of non-condensable gases. Also, a method was developed 
to evaluate thermal resistance of the condenser tube walls due 
to the axial flutes. 
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I . INTRODUCTION 



A design optimization procedure for vertical steam 
condensers is established here to develop condensers with 
minimum weight and volume . The performance comparison of 
enhanced steam condensers conducted in Reference-3 suggests 
the following methodology that provides the basis for this 
optimization scheme. 

The optimization procedinre is outlined as follows : 

1. Design for the following given condenser operating 
conditions: 

a. total steam condensed 

b. condenser operating pressure 

c. coolant inlet temperature 

d. coolant flow rate 

2. This design is subject to the following specified 
parameter limits: 

a. maximum condenser length (L ) 

^ ' max 

b. maximum pumping power 

c. coolant velocity range from a minimiM velocity 

(Vjjjin) to a maximiira velocity . 

3. For a given heat transfer surface minimum weight 
and volume will always be achieved for a maximum coolant flow 
velocity that results in; 

a. condenser length is equal to for pumping 

power less than P ,, or 

max 

b. pumping power equal to Pjjj^ for condenser length 

less than 

max 

Thus, the optimization routine proceeds as follows: 

a. select value for coolant velocity (V) 

b. calculate overall heat transfer coefficient (U) 
using the equations from Chapter II 
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c. calculate tube length and the rest of the 
condenser geometry to finally determine total condenser 
length (L ) 

d. check if pumping power (P) is less than 



and L = 
c 



L 



max 



e. if the above conditions for P and L are not 

1 ^ 

met, then either change V or enhancement of heat transfer 
surface 



5 . a. 



and P<P 
or Pa P 



max 

max 



if the calculated condenser length L_<L 

c IH3.A 

change to less enhanced surface until L^a 



b. if L and F<P increase the surface 

c max max 

enhancement until either L^= L 

c max 

c. if L >L__„ and P>P„„^ 

c max max 

either L a L and P<P or Pa P_„„ 3-nd 

^ max max max max 



OF 

* max 

then decrease V until 



Further details of this optimization procedure can be 
found in Chapter V-Computer Modeling, 

Vertical condenser tubes can provide great flexibility 
in this design optimization procedure due to the different 
possible combinations of internal and external enhancement 
configurations. Thus it is possible to "fine time" the 
overall heat transfer coefficient to obtain the proper con- 
denser length. 

The procedure that has been outlined here is included 
in a preliminary design program VERTCON-2 along with the pro 
posed vertical condenser arrangement shown in Pigure-1. 

Steam enters the condenser shell along a major portion 
of the tube length, steam lanes distribute the steam aroimd 
the tube bundle with radial steam inflow into the tube 
bundle towards the air cooler section and air removal duct. 
Steam condenses on the vertical tubes , and the condensate 
collects in the valley of the tube fluted surface flowing 
down the tube length. The tube support plates have a dual 
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FIGURE 1 

VERTICAL EMHAIMCED TUBE CONDENSER 
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purpose by also serving as condensate drain-off plates. 

The condensate is stripped off the tubes and collected on 
the tube support plates, and then it is removed to the hotwell 
section via downcomer drainage tubes. Improved condensate 
control is afforded by this condensate drainage and removal 
scheme. The tube support/condensate drain-off plates 
divide the condenser up into sections, where the full depth 
of tubes in each section see relatively fresh steam. The 
tubes will see fresh steam through the depth of the tube 
bundle since the tubes are not subjected to condensate 
indundation effects because of effective condensate removal 
and vertical orientation of the tubes. A continuous duct 
that runs along the total tube length is utilized for air 
removal. Figure-1 shows a double flow coolant water circuit 
with the inlet/outlet header located in the bottom of the 
condenser. The condenser hotwell is located below and also 
encases the inlet/outlet header. This hotwell location adds 
to the total height of the condenser. 
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II, HEAT TRANSFER CALCULATIONS 



Nomenclature 
A area (ft^) 

a amplitude of the flute (ft) 

Cp specific heat (BTU/lhm®p) 

D»d diameter (ft) 

e helical ridge height (ft) 

f friction factor [■(AP/p)(D/L)(2g^/v^)] 

G mass flux (Ibm/ft^hr) 

g^ gravitational constant (Ihm ft/lbf hr ) 

h heat transfer coefficient (BTU/hr ft^**?) 

h^g latent heat of vaporization (BTU/lbm) 

K thermal conductivity (BTU/hr ft^ ®?) 

L length ( ft ) 

1 lead of ridge (axial distance per 360®turn)(ft) 

m operand for friction factor equation 

Nu Nusselt number 

2 

P pressure (Ibf/ft ) or (in-hg-abs) 

p pitch of ridging or flute (ft) 

Pr prandtl number 

Q heat flow (BTU/hr) 



Re 

Re^ 

% 

^s 

St 

T 

AT 

^"^Im 

U 



operand for friction factor equation 
Reynolds number 

flooding Reynolds number for flute |^= ** 

thermal resistance of axial flute (hrft^ ®?/bTU) 
thermal resistance of scale (hrft^ ®p/bTU) 

thermal resistance of tube wall (hrft ^ ®f/bTU) 

S tanton number 
temperature (®F) 



mean overall temperature difference ( ®F) 
overall heat transfer coefficient (BTU/hrft^ ®F) 
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U*^ Dimensionless velocity [■u/u* ] * 

U* friction velocity j * 

V velocity (ft/hr) or (ft/sec) 

V/- flooding axial mass flow of condensate per flute 

^ ( Ibm/hr ) 

half-perimeter length of flute (ft) 

Subscripts 

b fluid at the bulk temperature ( *P) 

c condensate 

f fluid or flooding 

i inside, or inlet 

n nominal 

o outside, or outlet 

s increment, or section 

sat saturation 

w wall 



Superscripts 

+ dimensionless parameter 



Greek Symbols 

a 



7 

X 

M 

P 

<T 

Q 

T 



height of the condensate in the center of the 
flute (ft) 

fragment of -2.51n(2e/d£) 4 3*75] * 



dimensionless group s 4 p' 



g. _Lfo) 



•»« 



dynamic viscosity (Ibm/hr ft) 
density (Ibm/ft^) 
surface tension (Ibf/ft) 
non-dimensional group -f* 
apparent wall shear stress 



* Reference (6) 
** Reference (4) 
f Reference (1) 
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II. A ANALYSIS SUMMARY 



The rate of heat flow Q over an entire heat exchanger 
is related to the mean overall temperature difference 
and the total heat transfer area A hy the overall heat 
transfer coefficient U. 

Qs UA (II-l) 



From the analysis of externally ridged tubes it proves 
convenient to base the heat transfer coefficient U on the 
surface area of a smooth tube having an outside diameter 
^n diameter measured over the mid-height of 

the external flutes (see Figure- 2 for external flute 
geometry) . 

V — ^ (II-2a.) 

tD„L 

where: D„= + 2a (II-2b.) 



The overall U depends on the resistances in series ; 
between cooling water and the tube wall, within the tube 
wall, and between the tube wall and the working fluid. 
Fouling resistance on either side of the tube wall can 
be combined in one term R^. 

Tube wall resistance is considered to be the summation 
of an equivalent smooth tube resistance and flute resistance 
Tube wall resistance for a smooth (non-enhanced ) tube can 
be expressed as: 



V 1 w ln(Do/0^)D^,f 



(II-3) 



Where is the reference diameter on which the overall 

U is based on. For a fluted tube the overall U is based 
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on the nominal diameter D^. As will be discussed in the 
next section , the flute resistance can be written: 



R 



flute 



= lnf7.985(a)-7^^^(p) 



-.0242 



4irref 



2k 



(II-4) 



Therefore, the total tube wall resistance for the fluted 
tube can be written: 

- . 0242 

lnK/Dj^7.985(a)-'^^^^(|) -*-lil^n ( 11 - 5 ) 
2k 

Therefore, the concept of resistances in series yields: 




h. 



'Di 



Rs +Rw 



1_ 

h. 



(II-6) 



The overall heat transfer coefficient in Equation 
(II-6) does not consider non-condensable gas effects. 
Condensing in the presence of non-condensable gases will 
be taken up in Chapter V. 




FIGURE 2 

Cross section of Vertical Fluted Condenser Surface 




Z -coordinate into Page 




Application of Heat Flow Restriction 
Model to Axial Fluted Tubes 
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II. B TUBE WALL RESISTANCE 



One method to evaluate tube wall resistance for fluted 

tubes would be to use Equation-3 » where and equal 

the nominal diameter D^. 

n 

R s in(Dn^Dj)Dn (II-7) 

^ 2K 

This formulation is not sensitive to flute geometry 
such as changes in amplitude -to -pitch (a/p)ratio. In com- 
paring Equation-7 with finite element formulation differences 
of up to 30 percent in flute resistance are observed. This 
could amount to 10 percent difference when considering 
total tube wall resistance for tubes made with low conduct- 
ivity materials such as titanium with thin tube wall dim- 
ensions . 

What is proposed here is to model flute and condensate 
layer as a heat flow restriction, by describing two-dimensional 
heat flow through an elliptical boundary. (Figure- 3&) 

This serves as an accurate model for a heat flow restriction 
as shown in pigure-3b. 

The procedure is to apply this same model to describe 
the heat flow through the axial flutes (Figure-3c). 

Some aspects of curvilinear coordinate systems (7) 
will be reviewed in order to develop nomenclature to be 
used in later equations . 

Rectangular coordinates x,y,z, can be expressed in 
terms of new coordinates u^, Ug, u^ by the equations; 

(see Figure-^) 

xs x(u^, U 2 » u^) 
ya y(u^, u^) 

z= z(u^, u^) 
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The coordinates are chosen such that u^, u^, are 

orthogonal curvilinear coordinates. From exaunination of 
Figure-^ it can he seen that if a particle moves from 
point P in such a way that 1 I 2 and u^ are held constant 
and only u^ varies, then a curve in space is generated. 

Let r represent the position vector of a point P 
in space: 

— A A A 

rr xi 4 - yj ^ zk 

where 1 , and are \mit vectors in the x, y and z 
directions respectively. 

Then a tauigent vector to u^ curve at point P is given 

by: 

V- N- dS 

Ul* = -i 1 (II-8) 

^U]_ ^S3_ du^ 

arc length along the u^ curve 
^ is a unit vector 



thus : ^1^1 



where u^is the unit vector tangent to u^ in the direction 
of increasing arc length. 



thus : 



h, - ^^1 which is the length of U, 
du^^ ^ 



Similarly for the other coordinates : 

U2=h2U2 U3= h3U3 

hi= f£i is 1.2,3 dS^M h^du^ 



(H-9) 

i= 1.2,3 
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^^3 

- y 

FIGURE ^ 

Curvilinear Coordinate System 





; Steady One -Dimensional 
; Heat Transfer Through 
.* A Thickwalled Cylindrical 
; Shell. 
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The unit vectors u^, and form a right-handed 
orthogonal curvillinear coordinate system. The vectors 
UidUi* U2 U^du^ are mutually perpendicular vectors 

having as their length the arc-length differentials dS^, 
dS2 » and dS^. Therefore, an element of volume can be 
determined from the following vector product: 



dVs U^du^ X U2<iU2 
dVs h^h2h^du^du2du 


• U^du^ 

(II-IO) 

3 


Also, a vector element of surface area dA, on the surface 
u^r constant is given by the vector product: 


dA2« h2h^dU2du^ 


(Il-lla) 


and similarly: 




dA2= h^h^du^du^ 


(b) 


dA^r h^h2du^dU2 


(c) 



GENERAL PROBLEM 

Consider steady one -dimensional heat transfer through 
a thick-walled cylindrical shell a section of which is 
displayed in Figure- 5* 

First assume rate of heat transfer is constant at 
every section. 

Qs constant 

Q«“KA^ (11-12) 

ds 



Where S is defined in Figure-5, integrating Equation (11-12) 
between inner and outer surfaces with the assimption K is 
constant yields : 



Q= 



T1-T2 






= ^ (II-13) 



where R is heat transfer resistance. 
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thus : 



R= 



I = i /* 
Ky 



_^T = 
Q 



ds 



dA(s ) 






Equation (II-14) is only applicable when heat flux is con- 
stant across the heat transfer area. In many geometries 
the distribution of heat flux across the heat transfer area 
may not be uniform, but equation (II -14) can be expanded 
for variable heat flux ( 16 ) . 

Equation (II -14) can be rewritten: 






ds2ds^ 



(11-15) 



where : ds^= h^du^^ 

dS2* ^2^^ 
ds^= h^du^ 

thus : ^ = i dT 

dsj^ hj^ dUj^ 



and substituting into Equation (11-15): 



1 A 



thus : 



Equation (II-16) is applicable to variable heat flux dis- 
tribution over the area of heat transfer. But this equation 
is restricted to one-dimensional , steady state analysis, 
where temperature is a function of u^_ only. 

Equation (II-16) can be applied to the model of two 
dimensional heat flow through an elliptical boundary. 

Referring to Figure- 3 and Equation (II -16) the following 
expressions are obtained : 



R= AT 
Q 



(Ut ) 




(11-16) 



(Ui)i 
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thus : 



7^r constant temperature lines 

7,^jZ coordinates , where z coordinate is 
into page 



^i= 



hr 

au. 



where r r 



XI 



yj 



A 

zk 




(II -17) 



in this case: 

x= c coshn COS’/' 
ys c sinhn sin>/' 
Zr z 



from Equation (II-l?)* 

^ 1 = ^2 
h^« 1 



substituting into Equation (II-16) 
R= 1 / d7 



t= if _ 

r, 



If 

0 



I d^dz 



R= ’^l 
KxL 



and : 



(11-18) 
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x= c cosh’7 cosiA 
at <Az 0 

Xs t) r c cosh7 
and similarly: 
y- c sinh 7 sin^ 
at V2 
y = a = c sinhT? 

therefore ; 

V-^s sinh 



and from a fundamental identity; 



In 



c If 0 



+ 1 



(11-19) 



For axial fluted tubes Equations (11-18) and (11-19) 
will take the form; 



R= ^ 
K-^rL 



(II-20a. ) 



»? = In 






+1 



(b) 



where C is defined as the flute resistance coefficient. 

Equation (11-20) is in good agreement with actual physical 
situation ; 

If asO this implies 7s ln(l)zo, thus Rso 

1. Therefore a flute of zero amplitude has zero thermal 
resistance . 

2. As flute amplitude increases the flute resistance 
also increases. 

Equation (11-20) can be simplified as follows; 

1. Amplitude-to-pitch (a/p) ratios for fluted tubes 
are in the range of .l-to-.5* 

2. The flute resistance coefficient (C) will have an 
approximate value of C- .5 thus, 
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and eqtiations (11-20) simplifies to: 

R= lJ^[c(§)-n] (11-21) 

KttL 



To express thermal resistance in units of (hrft f/bTU) 
multiply R from Equation (11-21) by a reference area 
This reference area should be the same reference area 
that the overall heat transfer coefficient is based on. 
thus : 

Ra fref^f^ (hrft^ ®F/BTU) (11-22) 

Q 



In 






R^* ^J^ref (11-23) 

2K 



Values for flute resistance are calculated from a finite 
element formulation using the ADINAT (2) finite element pro- 
gram. The flute resistances obtained from ADINAT are sub- 
stituted into Equation (11-23)# and the values for flute 
resistance coefficient (C) are determined at different 
values of a/p, and for various tube dimensions . One result 
is that C is a function of both flute amplitude a, and 
a/p ratio. Thus C as a function of a and a/p was correlated 
using multiple non-linear regression ( 15 ) analysis. 

Then a final formulation for flute thermal resistance 
was established in the following forms 




(11-24) 
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^ , f 

where C= \pj 

^o’ ^1’ ^2 correlation coefficients resulting from the 
multiple non-linear regression. 

The convection heat transfer coefficient along the 
perimeter of the flute is required for input into ADINAT. 
The distribution of convection heat transfer coefficient 
was estimated by first calculating the condensate film 
thickness (5) along the perimeter using a procedure estab- 
lished by BARNES (^). The convection heat transfer at a 
particular point along the flute perimeter can be approx- 
imated by: 

h~ T (BTU/hrft^ V) 

where K is the thermal conductivity of the condensate, and 
Sis the condensate film thickness. 

As stated before, the finite element results, along 
with regression analysis were used to determine the value 
of C shown below : 

C= (11-25) 

b^- 7.985 
b^_r .7688 
b2=“ 1.0242 

Finally, equation (11-24) can be written; 

^ 2K 



The final results of this analysis is depicted in 
Figure (6). Curve-1 represents flute thermal resistance 
as a function of a/p ratio for constant amplitude flutes. 
In Curve -2 both flute amplitude and a/p ratio are varying, 
but flute pitch is held constant. 




I 
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Also , compared on this plot are the flute thermal resistance 
calculated from the finite element results , from Equation 
(II-7) based on nominal tube diameter, from Equation (11-26) 
based on heat flow restriction model. For ease of comparison 
thermal resistance have been normalized on the maximum 
range of thermal resistance considered, which is denoted 

?igure-l6 shows that the discrepancy in value between 
the finite element solution and result based on nominal 
tube diameter from Equation (II-7) increases with increasing 
flute aimplitude and a/p ratio. Both Curve-1 and Curve-2 
together cover a broad range of flute amplitude and a/p 
values to demonstrate the fit of Equation (11-26). 

From the results that were depicted in Figure-6 indicate 
that Equation (11-26) accurately models the flute thermal 
resistance. This equation is based on finite element sol- 
ution which must also be checked to ascertain whether it is 
a valid model. If the finite element model is developed 
methodically from a simple configuration where the results 
are known, to the more complicated model, then confidence 
can be placed on the results. In this case accurate spec- 
ification of convection heat transfer coefficients was 
important to the accuracy of the finite element solution. 
Also, a simple check can be made to determine if the finite 
element grid has the proper number of elements to model the 
problem. This is done by checking the heat flux approaching 
a boundary node from different elements. If the grid is 
sufficient there should be no jump in the value of total 
heat flux across the bo\mdary at a specified location. 

These general guidelines were followed in modeling tube 
wall resistance for axially fluted tubes. 
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Curve-1 - Constant a (as .002 ft.) 

Curve-2 - Constant p (ps .OO 625 ft.) 

■■ — Finite Element Solution 

Equation (11-26) Heat Flow Restriction Model 

Equation (II-7) Average Flute Resistance 

R*s 4.03 X 10"^ (hrft^®F/BTU) 



1 . 0 - 

.9 

.8 

.7 

Rf .6 

.4 

.3 

.2j 




I I I T" 

I .2 .3 .4 

a/p 

FIGURE 6 



Comparison of Flute Thermal Resistance Calculations 
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II. C WATERSIDE HEAT TRANSFER COEFFICIENTS 



1. Smooth Internal Tube 

For the smooth internal tube the McAdams correlation 
was used for determination of the heat transfer coefficient 
for cooling water. 




0.023 



GD ■ 


0.8 


‘mCP 


^ J 




k 



0.4 



(II-27a. ) 



or : 



h- 0.023 



0.4 



(II-27h. ) 



where ; 

(a) all fluid properties are evaluated at the bulk 
fluid temperature? 

(b) 2300 Re^^ lo"^ where; 

Rej^s Reynolds number based upon hydraulic diameter 

(c) 0.5 Pr.j^ 120 where; 

Pr^= Prandtl number based on bulk temperature 



2. Helical Internal Ridging 

Internal tube augmentation was also investigated, with 
the application of integral multiple helix ridging ( 13 ) as 
sho’ym in Figure- 7* 

Friction factor data was correlated by means of the 
following equation: 




1 

2.46 ln[r+(7/Re)“^] 



(11-28) 



For the roughness of the helical internal ridge both r and m 
are treated as variables, and are tied in with tube geometry. 
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FIGURE 7 

Internal Enhancement Geometry 
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Both in and r vary with the dimensionless parameter e/l, 
where e is the ridge height and 1 is the lead of the ridge. 
In reference (13) WITHERS has made a distinction based 
on the criterion p/d se O .36 in correlating the friction 
behavior to the internal geometry. It has been proposed 
that a shift in flow behavior occurs at p/d m 0 . 36 . 

For higher values of p/d a greater degree of swirling could 
occiar, compared with cascading of flow that takes place if 
p/d is less than O. 36 . 

The heat transfer correlation equation developed from 
data for tubes of various configurations , and solved for 
the Stanton number becomes: 



st= ^ ^ 

5 . 68(e/p)"^'^“T/pr [(e/d^)Re^/f78]°'^36^ y 

where: /r “[2-5 ln(2e/d^) -#>3.75] (II-29b.) 



(II-29a. ) 



Then, h^ becomes: 




K Re Pr St 
c 



(II- 30 ) 



These equations are applicable to Reynolds number 
range 10,000-120,000 and Prandtl number range ^-10. 

See Appendix-B for tube data, and friction factor 
characteristics of multiple -helix internal ridged tubes. 
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II. D CONDENSATE HEAT TRANSFER COEFFICIENT 



In reference-2 BARNES developed an equation for the 
average value of condensate heat transfer coefficient (h^) 
for an externally fluted tube of length L. In this 
formulation h depends upon physical properties of the 
condensing fluid geometric factors, and mass flow rate of 
condensate in the flutes. 



\ z .6027 




.9226 




.2307 
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The non-dimensional group Nu^n** defined in Reference-2 is 
a function of flute amplitude -to -pitch (a/p) ratio. It 
should be noted in Equation- 31 the parameter L is the 
length of the tube between condensate drain off plates. 

The parameter is the flooding axial mass flow of 
condensate per flute (LBM/hr) which is also a function of 
physical properties of the fluid, aind tube geometry. In 
Reference-?, PANCHEL and BELL defined the following flooding 
Reynolds number for condensate flow in the axial flutes : 

RSfS fwi ( 11-32) 



and an additional non-dimensional group: 



A 





( 11-33) 



where “o- 2a in the case where the flute is flooded. Also, 
the following correlation was developed for A^ based on the 
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half -perimeter length (Xj^) of the flute; 

36(a/p)exp(3.33 a/p) Re^- (11-3^) 



Substituting equations (32)and 33)iri'to equation (3^) and 
solving for yields : 



8 a^p exp [-3.33(a) 

i 9 M p 



(11-35) 



where is now based on full perimeter length of the flute. 

The non-dimensional group Nu^n** can also be defined 
by the following correlation; 

Nu.Q^ * h 4 f] b,(a/p)^ (11-36) 

° ° k-1 ^ 



where b^, bj^ are the correlation coefficients; 
b = 3.661218 
b^s 6.^526 

b2=“15.265 

b^= 16.1^5^3 

bj^=~6. 561166 

Equation (11-36) was obtained from multiple linear 
regression, using techniques of Reference-15. 
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II. E SELECTION 0? EXTERNAL ENHANCEMENT GSCMSTFY 



The heat transfer performance of various vertical tubes 
with different levels of external enhancement were examined. 
This was done to provide guidance in the selection of external 
enhancement geometry. The basic parameters that need to 
be specified are: 

1. Flute amplitude (a) 

2. Flute amplitude-to-pitch (a/p) ratio, 
and thus flute pitch (p) is also set. 

3. The ratio of condensation rate-to-flooding 

condensation rate . 

The following plots of h /h * versus a/p ratio (Figure-8) 

illustrate the relationship between these basic parameters. 

The variable h^ is the condensing heat transfer coefficient, 

which is normalized in the plots by the maximum value of 

h (denoted h *). The curves of Figure-8 denote lines of 

constant flute pitch (p) where curve-1 is for w/Wp= 0.1 

and curve-2 is for V//w^s O.25. it can be observed that 

^ . 
greater enhancement is achieved with larger values of a/p 

ratios, and smaller values of flute pitch. In Figure-9 

a line of constant flute amplitude is also plotted at w/w^=.l. 

Pigure-9 shows that increased enhancement is more readily 

achieved from increasing a/p ratio by reducing flute pitch. 

Increased steamside enhancement cannot be achieved 
without paying a cost. The cost in this case is increased 
tube wall resistance. In Figure-10 normalized flute resistance 
for constant amplitude flutes (curve-1, a= .0125”) and 
constant pitch flutes (curve-2, p® .025") are plotted along 
with heat transfer data for w/Wp- .1. Figure -11 is a 
similar plot except curve— 1 is for ar .01, and curve— 2 is 
for p* .050. Flute thermal resistance was calculated using 
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FIGURE 8 - constant Flute Pitch Curves 



a=.0125 in. 
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FIGURE 10 - Constant Amplitude (as.0125in.) and Constant 
Pitch (p= ,025in.) Curves 
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Equation (11-26). From Figures 10 and 11, the following 
conclusions can be inferred : 

1. V/hen comparing on the basis of flute pitch 

at the same a/p ratio, it can be seen that flute resistance 
is larger and heat transfer enhancement is less for 
increased values of flute pitch. 

2. v.’hen comparing on the basis of flute amplitude 
at the same value of p, it can be seen that flute resistance 

is larger and heat transfer enhancement is greater for increased 
values of flute amplitude. 

This indicates that the best tube geometry is small 
amplitude flutes with larger a/p ratios , and thus smaller 
values of flute pitch in comparison with other tubes. Thus 
the benefits of increased heat transfer enhancement are ob- 
tained with minimum increase in flute thermal resistance. 

As demonstrated by Pigure-8, greater heat transfer 
enhancement is achieved at smaller values of Vv/v/^, but there 

r 

is a lower limit in practical condenser design for 
Reduced values of condensate drainage plate spacing are 
required to achieve smaller values of Close drainage 

plate spacing will complicate condenser fabrication and 
maintenance. Therefore acceptable values of condensate 
drainage plate spacing should be used as criteria for specify- 
ing W/Wp in a practical condenser design. This criteria is 
used in the program VERTCON-2 to specify This effect 

of on condensate drainage plate spacing can also be 

observed in pigure-8, where the drainage plate spacing at 
Point A (wA’ s .25) is equal to 1.632 ft. while at Point B 

r 

(W/Wp= .1) the spacing is 0.622 ft. 
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III. CONDENSATION IN THE PRESENCE 0? NON-CONDSNSABLB C-ASSS 
Nomenclature 

2 2 

A Flow area (m ) , (ft ) 

a Constant defined in Equation (III-6) 

B Term defined by Equation (III-l4a) 

Cpf Specific heat of cooling water (j/kg^^K), (BTU/lbmV) 

D External diameter of tube (m), (ft) 

o 

Nominal diameter of axial fluted (m), (ft) tube 

D Coefficient of vapor diffusion in gas (kg/m s)/pa, 

^ (Lbm/ft s)/Lbf/ft2 

2 

G Vapor mass velocity approaching tube rov/ (kg/m s), 

(Lbm/ft^s) 

2 2 

g Gravitational acceleration (m/s ), (ft/s ) 

Specific enthalpy of vapor (j/kg), (BTU/Lb) 
h^ Specific latent heat (j/kg), (BTU/Lb) 

Specific enthalpy of vapor (j/kg), (BTU/Lb) 

L Length of condenser (m), (ft) 

N Number of tubes per row 

Partial pressure of air in bulk of mixture (P„), (Ib/ft'^) 

3 . p d 

P_ Pressure of mixture (P„ ) » (Ib/ft"^) 

Ul d 

P„ Partial pressure of vapor at condensate surface 

® (p^). av«2) 

P Partial pressure of vapor in bulk of mixture 

^ (P^), (lb/ft2) 

q Heat flux (W/m^), (BTU/hrft^) 

Temperature of coolant (*K), C*R,“P) 

Tj^jyj P/lean temperature of coolant in tube ('’k), (°R,'’f) 

T„ Temperature of condensate surface (‘’k), (°R,*’f) 

w 

T^ Temperature of vapor (*’k), (“r.^F) 

T,, Temperature of outer tube wall (°K), (*’r,^F) 

w 

U Heat transfer coefficient for convection and tube « 

resistances as defined by Equation (II-6) (BTU/hrft^ ®F) 

U Overall heat transfer coefficient between vapor and 

coolant (W/m2*K), (BTU/hrft^ '’ f) 

Specific volume of vapor (m^/kg) , (ftVlb) 
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Flowrate of air (kg/s), (Ib/s) 

3 . 

W Rate of condensation per row (kg/s), (Ib/sec) 

c 

Flowrate of cooling fluid (kg/s), (Ib/s) 

Flowrate of vapor (kg/s), (Ib/s) 

w 

s Condenser steam flow (Kg/hr), (Lb/hr) 

Subscripts 

N Subscript indicating nominal tube diameter 

V Subscript indicating vapor 

1 Subscript indicating conditions at coolant inlet 

2 Subscript indicating conditions at coolant outlet 



Greek Symbols 



^w 

*^Wo 



volume ratio of air to mixture 

Weight ratio of air to mixture 

Weight ratio of air to mixture at outlet of air 
cooler 

Absolute viscosity of air (n/m s) 

Absolute viscosity of mixture (n/m s) 

Absolute viscosity of vapor (n/m s) 

Term defined by Equation (III- ) 

Coefficient of mass transfer (kg/m s)/p 



Dimensionless Groups 

NUjj Nusselt number defined by Equation (III -6) 

8®jjj Reynolds number of vapor gas mixture 



HEAT AND mSS TRANSFER EQUATIONS 



A significant decrease in condensing heat transfer 
coefficient (h ) can result from the presence of very small 
amounts of non-condensahle gases. The motion of the vapor 
with entrained gases moving towards a heat transfer surface, 
results in an accumulation of non-condensable gases. A 
vapor partial pressure gradient is required for the vapor 
to be able to diffuse through the gas blanket. This results 
in a liquid-vapor interface temperature below the temperature 
of the main vapor gas mixture. The heat and mass transfer 
equations developed by CHISHOLI^ (6) were applied in the 
VERTCON-2 condenser design program in predicting the con- 
densation of steam in the presence of non-condensable gases. 
Diffusion coefficients obtained from experiments were used 
in this analysis by CHISHOLM (6). 

The motion of the vapor towards the condenser surface 
is described by: 

The rate of heat transfer per unit area is given by: 
q= (III-2) 

and, if its assumed that only laten heat is removed, them 
G^h^g (III -3) 

Prom Clapeyron's approximate equations assiuning that 
P^-Pg is small, then: 



-40- 



By combining equations III-l, III-2, and III-4 results 
in the governing heat transfer equation within the vapor- 
gas mixture shown below: 

h - 2 



The coefficient of mass transfer is evaluated from: 



NUj^r (III-6) 



where* 



NUp^r ^*^^0 
°p 



(ni-7) 



Rem= -yf.a. 

A 



(III-8) 



'm 



1 = 



P -P 
V s 



m 



(III-9) 



For R®m'>350 

a= 0.52 for first tube row 
= 0.67 for second tube row, and 
s 0.82 for third and later tube rows 

For RSj<350 

ar 0.52 for all rows 



The absolute viscosity of the mixture is determined from 
'‘m= ^ l-6.£^»‘a]/(l+0.6l£^) (III-IO) 
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The coefficient of steam in air is obtained from? 



D^= 2.93x10 
P p 

m 



-9 



Zv 

L273 



0.8 



(III-ll) 

(kg/m s)/Pa 

The ratio of air volvime to total mixture volume is obtained 
from: 

< HI-12) 

and the partial pressures are evaluated from the equation: 



Pm € 

— - 1 + 0.622 W 






Combining equations III-6, III-7* and III-9 yields 



qr B(T^-Tg) 



2/3 



(III -14a. ) 



where. B= ^ ^ p___l/3 



T V 

. V V. 



2/3 



h^g^/^ III-l4b.) 



The heat transfer coefficient between the condensate 
and the cooling water is evaluated from: 

q= U(Tg-T^) (III-15) 

where U is the overall heat transfer coefficient defined in 
equation (II-^). 

Combining equations III-14 and III-15 yields: 



rnl3/2 r.-i 



£ 

B 



= (III-16) 



qr (T^-T^)/ 



B 



+ i 

yz u 



or : 



(III-17) 



— 



For given values of B» and U and temperatures and 
this equation can be solved for q. Equation (III-l?) 
will converge for any positive initial value of q. For 
the first tube row and first iteration assume q is zero 
and the coolant temperature is the inlet value, and for 
the next row and first iteration assume q and mean temp- 
erature of the coolant from the previous row. 

In evaluating U the condensate properties are first 
evaluated at the vapor temperature. After obtaining q on 
this basis , the overall heat-transfer coefficient between 
vapor and coolant is calculated from: 

UqS q/CTy-T^) (III-18) 



Using the overall heat transfer coefficient and log 
mean temperature difference for the coolant, Equation 
(III-18) can be written: 









(III -19) 



substituting: Qs (III-20) 



results in: 

T *T 

= expC-U^TrOj^LN/W^Cp. (III-21) 

V f 2 

and the mean temperature of the coolant is evaluated fromi 



''iM = (111-22) 

For successive iterations is used as the cooling 
fluid temperature (T^) in Equation (III-l?). At this point 
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the outer tube wall temperature and surface temperature of 
the condensate can be calculated, so that a more accurate 
value of U can be determined. This procedure is repeated 
with more accurate values of U, B, and until successive 
values of q are within specified tolerances. 

Then, the total condensation for a tube row is cal- 
culated from: 



q^DjjLN/h^g (III-23) 

Standard design methods use either the BEAMA (1?) 
code in the UK or the HEI (18) code in the US. Both codes 
recommend that the air ejector should be sized to reduce 
the vapor gas mixture leaving the condenser to 7.5 below 
the saturation temperature corresponding to condenser 
pressure. As an example a vacuum of 29 in.Hg corresponds 
to a ratio of non-condensable gas to mixture flowrate 
to the ejector of: 

Figure- 12 can be used to determine the flowrate of non- 
condensable gases (W^) to be used in the sizing of the air 
ejector. In this study the BEAMA code was correlated using 
multiple linear regression (15)» resulting in: 

W. s 10.0 + . 182 3x10'’^ -.975x10"^% ^ (III-24) 

w S 



ORVAIR DUTY - Ib/h 
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FIGURE 12 



Recommended Air Extraction Capacity 



IV. PRSSSUP.S DROP CALCULATIONS 



Nomenclature 



A 

B 

C 

D, 






e 
E 

S 
L 

Ne 

P 
P 

^t 

aP 
aP 

en 

AP2 .AP3 

aPj, 



bm 

dm 



Minimum flow area between outer tubes of a tube 
bank (ft^) 

Constant defined by Equation (IV- 1 ) 

Constant defined by Equation (IV- 1 ) 

Conversion constant (Ib/ft to in.Hg) 

Equivalent diameter of flow area (the ratio of 
flow area to wetted perimeter multplied by four) 
(ft) 

Inside tube diameter (ft) 

Outside tube diameter (ft) 

Internal tube helical ridge height (ft) 
Dimensionless relative roughness = 

Gravitational acceleration (ft/sec^) 

Perimeter of tube bundle (ft) 

Condenser tube length (ft) 

Plow length along tube bundle perimeter defined 
in Figure- 13 (ft) 

Equivalent number of tube rows 

Pitch of internal tube helical ridging (ft) 

Tube pitch (ft) 

Total pressure (in.Hg) 

Specified condenser pressure (static) (in.Hg) 
Pressure loss through tube bank (in.Hg) 

Mean penetration pressure loss (in.Hg) 

Mean total pressure loss around tube bundle (in.Hg) 

Entrance pressure loss (in.Hg) 

Pressure loss along flow length L2 . L^, 
respectively (in.Hg) 




Waterbox inlet loss (ft. of water) 
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aP,^ 0 Waterbox outlet loss (ft. of water) 

APm_ Tube end loss (combined inlet and outlet), 
(ft. of water) 

plow velocity of cooling water (ft/sec) 

V plow velocity of steam (ft/sec) 

s 



Greek Symbols 

£. Tube surface roughness (ft) 

p Steam density (Ib/ft^) 

Steam viscosity (Ib/hr-ft) 



Dimensionless Groups 
f Priction factor 

f^ Priction factor defined by Equation (IV-I9) 

fgj^ Priction factor for internally enhanced tubes 

Re Reynolds number 






IV. A TUBESIDE PRESSURE DROP 

The total tubeside pressure drop is calculated from 
the summation of the following pressure drop components ; 

1. pressiare loss in condenser tubes 

2. tube end pressure loss (inlet and outlet) 

3 . waterbox pressure loss (inlet and outlet) 

The pressure loss for smooth internal vertical tubes 
is evaluated using a friction factor relationship developed 
by CHURCHILL (19) as shown in Equation (IV-1) which is 
applicable for any flow regime as a function of Reynolds 
number and relative roughness of the tube. 

f- 8.0 [(8/Re)^^+ 1.0/(A+B)^/'^]°-°®^^ (IV-1) 

where : 

A= [ 2.457 ln(1.0/(7/Pe)°®+0,27E)] 

(37,530/Re )^^ 

- surface roughness 

= 5xl0“^ft. recommended value (20) for smooth tubing 

2= dimensionless relative roughness where 

is the inside tube diameter (ft.) 

Pressure loss for multiple-helix internal ridged tubes 
is calculated using friction factor (f ) correlations 
evaluated from Equation (II -7) resulting from manufacturer's 
test data obtained from Reference-13. The friction factor 
data was correlated with f^^ as a function of Reynolds 
number, and the m & r operands with the following result; 




( 

I 
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2 




1 



(IV-2) 



2.46 In [r4-(7/Re)"‘. 



Tube properties and friction factor characteristics for 
multiple-helix internal ridged tubes can be found in 
Appendix -B . 

Tube end and waterbox pressure losses for single and 
double pass condensers are calculated using pressure loss 
correlations evaluated from multiple non-linear regression 
using HEI standards (18). Pressure drop data from the HEI 
code can be found in Appendix-B. The resulting correlations 
are shovm here : 

Single Pass condensers 
Waterbox Inlet Loss : 



P^I= .0l422v/-^^ 



(IV-3) 



Waterbox outlet Loss: 



p^0= .00237 v/-33 



(IV-4) 



Tube End Loss (Combined inlet & Outlet): 



PtE= -01787/ 



(iv-5) 



Double Pass Condensers 
Waterbox Inlet Loss : 



P^^r .01422V^^*^^ 



(IV-6) 



Waterbox Outlet Loss: 




(IV-7) 



Tube End Loss (Combined Inlet & Outlet): 



p^z .0356V / 



(IV-8) 



IV.B SHELLS IDS PRESSURE DROP 



Vapor condensing temperature is a function of vapor 
pressure f which follows a specific saturation temperature- 
pressure relationship. With the pressure drops encountered 
in the steam flow from the steam inlet to the air cooler 
section through the tube bundle, there is a resulting re- 
duction in condensing temperature in the direction of flow. 
This temperature reduction will affect the rate of heat 
transfer in different sections of the condenser. 

In Reference-6 HARRINGTON has noted that the condenser 
performance standards of the HEI (18) are based on static 
pressiure, but that the condenser tube more closely senses 
a temperature related to the total or stagnation pressure. 
The method for evaluating the steam pressure distribution 
established by HARRINGTON along with the assumption that 
condenser performance is based on stagnation pressure was 
used in this study. Also, this can be considered as a one- 
dimensional design method where the vapor flow distribution 
is represented by a single-flow path. The flow path is 
considered as cross flow, radially inward from the tube 
bundle perimeter to the air cooler section of the condenser. 

The pressure loss incurred by the steaun flow can be 
considered as two components: 

(1) steaun distribution loss for steam flow from the 
inlet to around the tube bundle 

(2) steaun penetration loss for steam flow through the 
tube bundle from the tube bundle perimeter to the air cooler 
section of the condenser. 
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steam Inlet 




FIGURE 13 

Steam Distribution Pressure Zones 
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Steam Distribution Loss 

Refer to Figxure-13for the following discussion. The 
total or stagnation pressure is calculated for pressure 
zone-1 and is evaluated from: 

Pt = ^ (IV-9) 

2g 

The steam velocity required in Equation (IV-9) and for 
the entrance loss calculation is assumed to correspond to 
the recommended maximum main steajn lane entrance velocities 
as shown in Table- 1 . 



Table- 1 



Recommended steam Design Velocities 



Condenser Design Recommended Maximum Main Steam 

Pressure (in.Hg) Lane Entrance Velocity, fps 

1 500 



2 



400 



3 

4 

5 



300 

250 

200 



The entrance flow area will have to be sized accordingly 
with respect to this entrance flow velocity. All pressure 
losses around the tube bimdle perimeter are based on the 
total pressure in zone-1. 

The next step is to calculate the entrance pressure 
loss (L^) from pressure zone-1 to the main steam inlet lane 
as follows ; 

2 

APen= Q ♦ Q5 (IV-10) 

2g 

The next step is to calculate pressure loss due to 
friction along L 2 . Calculate the steam velocity from the 
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flow area in this zone. It assumed the net steam flow is 
reduced proportionately along the flow length L2» L^, 
and L^. 

Reynolds number can then be evaluated from: 



Re= 



/^Vs 



(IV-11) 



where is the equivalent diameter of flow area, which 
is the ratio of flow area to wetted perimeter multiplied 
by four. 

Finally the pressure loss due to friction along L2 
is evaluated using the Tanning equation as follows : 




Dg2g 



(IV-12) 



The friction factor f is obtained from Figure- 14 , or 
by evaluating the following correlation obtained by using 
multiple non-linear regression analysis (15)» for an assumed 
roughness value. 

f= .0436 Re"'^®^ (IV-13) 

This procedure is repeated to determine the frictional 
pressure loss for and This calculation scheme can 

be conducted with a greater number of smaller length segments 
to obtain a more accurate pressure distribution around the 
tube bundle. 

The mean distribution pressure loss around the tube 
bundle can then be calculated as follows : 

Ll -fLg +Lj +Li^ 

This mean total pressure identifies the condensing conditions 
at the tube bundle perimeter. 
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steam Penetration Loss in Tube Bank 

Calculation of Pressure loss through a tube bank is 
accomplished by replacing the actual tube bank with an 
equivalent rectangular bank that has the same perimeter 
and number of tubes. The equivalent number of rows in the 
equivalent tube bank depth is given by: 



N 



e 



NP 

L 



(IV-15) 



Next calculate the velocity through the minimum flow 
area between the tubes. The flow area is obtained from: 



^f= (p-Dq^^t (IV-16) 



The steam velocity entering the tube bank can be cal- 
culated using this value for flow area (A). Then Reynolds 
number is evaluated from: 



Re= 



(Iv-17) 



The terms (/»&m) in Equation (IV-l?) are calculated at 
the temperature corresponding to the meain pressure at the 
tube bank entrance. 

The tube bank pressure loss equation, which is mod- 
ified for uniformly decreasing mass flow (21) 

2 

4f^(0.5 4 Ne ) C ^s (IV-18) 

3 2g 

The friction factor f^ is obtained from Figure-15 , or 
by evaluating the following correlation obtained by applying 
multiple non-linear regression analysis (15). 
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f^r 65. (lO^Re^lOOO) (IV-19a.) 

f^= .59lRe“-^^^(P/DQ)“^*^^^ (lOOO^Re^l.OxlO^) (IV-19b) 

HARRINGTON (6) has suggested that the following equation 
be used to determine the mean tube bank pressure loss: 

2 

4f^(0.5+__^)_£|Is_ (IV-20) 

The mean condensing pressure in the tube bank is then 
evaluated froms 

Pbm= Pt-"Pdm -"Pbm 

All of the heat transfer calculations are based on the 
mean condensing pressure (P-^^j^) “^be tube bank. 



FRICTION FACTOR, 



-55- 




FIGURE 

Friction Factor for Steam Distribution 

Losses 




FIGURE 15 

Friction Factor for Steam Penetration 
(Tube Bank) Losses 
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V. COfCPUTSR MODBLINC- 



The basic procedure for heat enchanger sizing that was 
employed in the VERTCON-1 program (3) is used in the Vertical 
Condenser Sizing Program No. 2 (VERTCON-2). In addition the 
VERTCON-2 program considers condensation in the presence of 
non-condensable gases, and the program optimizes the con- 
denser design for minimum weight and volume. Other features 
of the VERTCON-2 program are as follows: 

1. The option of smooth internal tubes or internally 
enhanced tubes. 

2. Design of single or double pass condensers. 

3. Option of submarine, or surface ship condenser 

design. 

4. Two basic configiarations are available for 
double pass condensers: 

(a) Conventional return header design 

(b) "U-Tube” construction 

5. Diagnostic messages used to warn when con- 
denser drainage plate spacing is beyond specifications based 
on heat transfer performance and mechanical design. 

6. Computation of recommended values for condenser 
drainage plate spacing and maximum main steam lane entrance 
velocity when they are not specified. 

Program Discription 

The program VERTCON-2 is described here , and diagram- 
matically in the simplified program flow chart (Figure-l6). 

The program is based on the FORTRAN-77 standard, and is 
listed in Appendix-A. Also included in Appendix -A are 
sample program input and output files. 

The program input consists of input data files, and 
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interactive input entered at the terminal. The input data 
file consists of tube geometry information for the various 
candidate tubes to be considered, and specification of 
condenser geometry, materials and operating conditions. 

The design features are selected interactively at the ter- 
minal. See Appendix-A for further details on data input. 

The program is set up for sizing just one condenser 
or a series of condensers. After inputing data the pro- 
gram prompts at the terminal for selection of design features. 
If condenser drainage plate spacing and maximum main steam 
lane entrance velocity are not specified the program cal- 
culates recommended values (Branch-A and Iteration-A, 
Figure-l6). The determination of these values are based 
on heat transfer performance and mechanical design to provide 
proper tube support. If these values are specified then 
the program checks if condenser drainage plate spacing 
is beyond recommended values (Branch-B , ?igure-l6 ) . If 
the values specified exceed those dictated by proper design, 
then diagnostic warning messages are printed at the terminal 
and in the program output. 

After input of interactive data at the terminal, the 
program conducts the selection of flow velocity and tube 
geometry to optimize the design as follows; 

1. For the first candidate tube geometry the 
condenser is sized at the specified minimum flow velocity 

2. (a) If condenser length (L) is less than specified 

maximiim condenser length ^ and pumping power is less 

than specified maximum pumping power then Step 1 

is repeated for a larger value of flow velocity. 

2.(b) This process is repeated until Ls L or 

P= P 



m.ax 



max 



3. If at V„. ;L then this candidate tube is 

min max 



rejected because of insufficient enhancement. 
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4. steps 1 through 3 are repeated for all can- 
didate tubes . 

The condenser sizing in Steps 1 through 4 is evaluated 
as follows : 

1. In iteration B and C (Figure-36) the total 
effective tube length for the condenser is calculated by 
adding N tube sections of length equal to the drainage plate 
spacing until a proper heat balance is obtained. Iteration B 
adds up N-1 of these sections of equal length. Iteration C 
then calculates the final section length, which may have 

a length less than the specified drainage plate spacing. 

2. After the total effective tube length is de- 
termined the condenser geometry, weight, volume, and system 
power are calculated for each condenser. 

It should be noted that the influence of non-condensable 
gases has not been considered yet, and a multiplicative 
factor (Fj^q) is simply used to account for non -condensable 
gas effects. 

At this point the minimum condenser sizing for each 
tube configuration are compared to select the condenser 
configuration with overall minimum weight and voliune. 

If the option is selected, the program will now consider 
the effects of non-condensable gases . Using the tube bundle 
and shell geometry of the selected condenser configuration, 
the mean pressure around the tube bundle is determined. 

With the conditions identified at the tube bundle perimeter 
the program steps through the condenser tube bundle row by 
row to calculate the effective tube length at local condit ions 
including non-condensable gas effects. 

If the condenser length calculated from the consideration 
of non-condensable gases differs from the original calculated 
value, then a new value for the multiplicative factor 
is selected and the complete process including tube and flow 
velocity selection is repeated until these length values 
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are equal. By conducting the calculations in this manner 
the non-condensable gas effects will be evaluated a minimum 
number of times, thus significantly reducing computer 
solution time. There is a limit to the number (Nj) of 
these iterations that will be completed by the program. 

If after iterations there has not been proper convergence, 
then warning messages are displayed on the terminal screen 
and printed in the program output. This indicates that one 
of the specified condenser parameters is out of limits for 
the particular condenser design. Non-convergence would 
most likely be due to factors that influence steam-side 
heat transfer performance. One culprit could possibly be 
too small a value for tube pitch-to-diameter ratio, that 
would cause excessive steam side pressure drop. 

The tabulated prograun output displays for each condenser 
the following information: 

(1) All input data and design features selected. 

(2) If applicable, diagnostic warning messages. 

(3) Detailed weight breakdown of condenser 
components . 

(4) Total condenser volume, center of gravity, 
and total weight in dry and wet conditions , and system 
pumping power. 







I 
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APPBNDIX A 



DATA INPUT 

See Table-3 for a sample input data file 
INPUT DATA FILE 

1. Assigned Logical Operator - 20 to specify input 
device. See lines 47-52 of progress listing for 
further details . 

2. Free Format 

3. Variable names starting with I, J, K, L, M, orN 
are integer variables and others are real variables , 
except for specific cases as noted. 

LINE 1 

Enter: SUB or SUR 

SUB - Submarine Condenser Design 
SUR - Surface Ship Condenser Design 

LINE 2 

Enter: NUMBT ,LMAX , PMAX, VT/IIN, VMAX 

NUMBT- Number of Tube configurations to be tested 
LMAX - Maximum condenser length (ft) 

PMAX - Maximiim pumping power (hp) 

VMIN - Minimum coolant flow velocity (ft/sec) 

VMAX - Maximum coolant flow velocity (ft/sec) 

Note t LMAX is a real variable 

LINS 3 

Enter: TCI, TSAT, QSTM, STMLD, PSAT PD 

TCI - Inlet Coolant Temperature ( F) 

TSAT - Steam Inlet Saturation Temperature ( F) 

QSTM - Quality of Inlet Steam to Condenser 
STMLD - Steam Condenser (Ib/hr) 
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TABLE-2: SAMPLE INPUT DATA FILE 
. .4 TUBE TEST CASE. . 



SUB 

4.24.34.275. . 8. .15. 

66. 1,1 43 . 89, . 85, 257 000. ,6. 5, 1.35 

2. . 6. .. 62 5, .5 62 5, .5 62 5, 2. 64, 488. 

558.. 488..488..488..558..282..1.56.7,.56,488. 
1 , . 555 , .62 5 , .01 , .245,0. ,0. 

48 . , 0 . , 0 . , 9 . 6 , .00033 

2. . 55 5, .62 5, .01, .245, .02 04, .191 
48. , .59 , .00197 , 9.6 , .00033 

2. . 555. . 62 5, .01, .245, .012 5, .47 50 
48. , . 762 ,0.0, 9. 6, .00033 

2. . 555. .62 5. .01. .245. .02 4. .0949 
48. , . 58, .0075,9.6, .00033 
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PSAT - Condenser Operating Pressure (in.Hg) Absolute 
PD - Tube pitch-to-diameter ratio 

LINE 4 

Enter; TIS, TOTS, TSHL, TTSP , THW , THDR , ITSD 

TIS - Thickness of internal tube sheet (in.) 

TOTS - Thickness Outer Tube Sheet (in) 

TSHL - Thickness of condenser shell (in.) 

TTSP - Thickness tube support plate (in.) 

THW - Thickness of hot well plate (in.) 

THDR - Thickness of Headers (in.) 

ITSD - Internal tube sheet density (Ib/pt^) (real variable) 
Note ; ITSD is a real variable. 

LINE S 

Enter; OTSD, SHLD, TSPD, HV/D, HDRD, TBD, ETUBE, TCOVER , COVERD 
OTSD - Outer tube sheet density (Ib/ft^) 

SHLD - Shell density (Ib/ft^) 

TSPD - Tube support plate density (Ib/ft^) 

HWD - Hotwell plate density (Ib/ft^) 

HDRD - Header density (Ib/ft^) 

TBD - Tube bundle density (Ib/ft^) 

ETUBE - Modulus of elasticity of tube material (Ib/in ) 
TCOVER - Thickness of U-tube header cover (in.) 

COVERD - Density of U-tube header cover (in.) 

LINE 6 

Enter; NTYPE, DI , DW, AE, AEPE, SI, PI 

NTYPE - Tube type; a. for smooth internal tube 

b. for doubly enhanced tubes 

DI - Internal tube diameter (in.) 

DW - Tube Wall diameter (in.) 

AE - External flute amplitude (in.) 

AEPE - Amplitude-to-pitch ratio of external ridging 

SI - Internal helix ridge height (in.) 
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PI - Pitch of Internal Helix (in.) 

Note ; For smooth internal tube (NTYPE-1) enter 00 for El, 
PI. Values of El and PI are obtained from tube 
data in Appendix-B. 

LINE 7 

Enter: FLUTE, MI, RI , KW , RSCALE 

FLUTE - Number of flutes per tube 

MI - Operand for friction factor, helical internal 
ridging (real variable) 

RI - Operand for Friction Factor, Helical internal 

ridging 

KW - Conductivity of tube wall (BTU/hr-ft- F) 
RSCALE - Fouling resistance (hr-ft- f/BTU) 

Note ; MI is a real variable. 

To consider more tube configurations, repeat lines 6 and 7 
for each additional tube to be included in the analysis. 



INTERACTIVE DATA INPUT 

Enter data at the terminal as follows : 

V/hen the program commences, the following message 
appears on the terminal screen: 



** WELCOME TO PROGRAM: VERTCON, ^^RSION-2 

ENTER PROGRAM RUN NUMBER 



Enter Run number, which can be an integer from O-to-999* 
Next at the terminal appears : 



ENTER THE NUTffiER OF CONDENSERS TO BE SIZED 
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Enter the numher of condensers to be sized, which is also 
an integer from l-to-999* 

Next at the terminal aooears ; 



ENTER "1" FOR SINGLE PASS, OR "2” FOR DOUBLE PASS CONDENSER 



Enter the proper integer mimber 1 or 2 for single or double 
pass condenser. 

Next at the terminal appears : 



Tt'JO BASIC CONFIGURATIONS ARE AVAILABLE FOR DOUBLE 
PASS CONDENSERS: 

(1) CONVENTIONAL RETURN HEADER DESIGN 

(2) "U-TUBE" CONSTRUCTION 

ENTER 1 OR 2 FOR CONFIGURATION SELECTION 



Enter the proper integer number 1 or 2 for configuration 
selection. 

Next at the terminal appears ; 



ENTER: (1) COOLANT FLOW RATE (GPM) 



Enter as real number with units of gallons per minute. 
Next at the terminal appears ; 



DO YOU WANT TO INCLUDE IN THIS ANALYSIS THE EFFECTS 
0? CONDENSATION IN THE PRESENCE OP NON-CONDENSABLE 
GASES YES OR NO 



Enter YES or Y, NO or N for proper choice. 
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Next at the terminal appears ; 



DO YOU WANT TO SPECIFY MAXIMUM MAIN STEAM LANE 
ENTRANCE VELOCITY (FT/SEC)? 

YES OR NO 



Enter YES or Y, NO or N for proper choice. 

If yes is entered, Next at the terminal appears ; 



ENTER MAXIMUM MAIN STEAM LANE ENTRANCE VELOCITY (FT/SEC) 

Enter this as a real number. 

If no is entered. Next at the terminal appears ; 



PROGRAM WILL SELECT RECOMMENDED MAXIMUM MAIN 
STEAM LANE ENTRANCE ^/ELOCITY USING CONDENSER 
OPERATING PRESSURE AS SELECTION CRITERIA. 



Next at the terminal aouears t 



DO YOU WANT TO SPECIFY CONDENSATE DRAINAGE 
PLATE SPACING (FT.)? 

YES OR NO 



Enter YES or Y, NO or N for proper choice. 

If yes is entered. Next at the terminal appears ; 
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enter CONDENSATE DRAINAGE PLATS SPACING (FT.) 



Enter this as a real number. 

If no was entered, Next at the terminal anoears ; 



PROGRAM WILL SELECT RECOMMENDED CONDENSATE 
DRAINAGE PLATE SPACING 



Note ; At this point, entry of Interactive Data is completed 
for the first condenser. The program now completes 
the sizing for the first condenser and will prompt 
the terminal for further interactive data for any 
additional condensers. 



TABLE-3 : SAMPLE OUTPUT DATA FILE 
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MATERIAL THICKNESS WEIGHT 

DENSITY(LB/FT**3 ) (IN.) (LB.) 

OUTER TUBE SHEET 558.0 6.0000 13879. 
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WRITE(KSCR , 108) 

F0RMAT(1X, 'ENTER THE NUMBER OF CONDENSERS TO BE SIZED 
READ(KSCR,*) NCOND 
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WRITE(KSCR , 122 ) 

122 FORMAT( IX ENTER "1" FOR SINGLE PASS, OR "2" FOR DOUBLE PASS CON 
IDENSER' ) 

READ(KSCR,*) NPASS 
IF(NPASS .LE. 1) THEN 



WRITE(KOUT ,124) 

F0RMAT(1X,'. . SINGLE PASS CONDENSER 

ELSE 

END IF 
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IF((GAS .EQ. NO) .OR. (GAS . EQ . NOE)) THEN 

NON= 2 

ELSE 

NON= 1 

END IF 
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KSAT= CONDUCTIVITY OF CONDENSATE ( BTU/ HR-FT-DEG . F ) 

MU= ABSOLUTE VISCOSITY OF CONDENSATE (LBM/HR-FT) 

SIGMA= SURFACE TENSION OF CONDENSATE (LBF/FT) 

SPVOL= SPECIFIC VOLUME OF STEAM AT SAT. CONDITIONS (FT**3/LBM) 



-79- 



U 

P-( 

w 

<3 

•K 

cn 

cn 

• 

m 



CO 



On 

X 

w 



w 

CLI 

w 

< 

w 

<3 



hJ 



^ W 
H 

CO ^ 
•JC hJ 

w 

<I »3 
•K CM 

eg w 
H 



4< 

O 

C/3 

2 

W 



< 

oc 

o 



CNI 

•K 

•jc 

W 

CM 

W 

<3 



uo 

vC 

eg 

• 

L/O <f 

I w 

CM 
W W 
CM <3 
W 
< 

•K vO 



vC 



























/-s 






Q 




eg 


SO 






















/-s 


/-s 




o 




eg 




Cxi 


uo 


uo 








H 




eg 


CO 


<r 


m 


vO 




00 


o> 








•X 






• 






eg 


oa 




























U 


• 


\D 








S 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


30 


H 


vO 




w 






33 






















15^ 




S 


< 




1 


Cxi 






2 












'w' 


V— ✓ 


>w/ 


'w' 


'w' 


'w' 


'w' 




CO 


+ 












<3 


<3 


<3 


<3 


<3 


<1 


<3 


<3 


<3 


< 


<3 


<3 




2 




CO 


II 




II 




H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 




w 


00 


■K 


/-s 




CM 




<3 


<C 


<3 


<3 


<3 


<3 


<3 


< 


<3 


<3 


< 


< 


uo 


Q 




•K 


CO 




o 




C3 


Q 


Q 


C3 


Q 


Q 


C3 


Q 


Q 


Q 


Q 


Q 


• 


2 


eg 


CiJ 




m 


o 


II 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


H 


m 


O 


*— 4 


CM 




eg 


hJ 


H 




























o 


vO) 


W 


H 


o 


2 


MC 
























II 


00 




vO 


< 




• 






II 














II 








W 




o 


• 








O 


m 


W 








II 






W 








hJ 




2 


CO 


•K 


<3 




eg 


eg 


CM 








W 






H 








<3 




M 








II 






>-• 


II 


II 


II 


CM 


II 


II 


33 


II 


II 


II 


O 


II 


Q 


II 




<3 


fX4 


O 


o 


H 


M 


::2 


W 


W 


M 


M 




M 


M 


32 


CO 


1X4 


O 


w 




Q 


S 


Q 


Q 


2 


C3 


Q 


<3 


<3 


w 


CM 




S 


CtC 






S 


O 


|X4 




H 






H 

2 

ClI 



w 

o 

o 

PiC 

w 

C/3 

2 

<3 

cc3 

H 

H 

<3 

W 

33 

O 

2 

M 

CO 

2 

W 

C3 

2 

O 



Ct3 

O 

Q 

2 

< 

CM 

O 



II 

[X4 

:2 



II 

w 



o o o o o 



o 



o 



o 



o o 



o o 



IF(AE .LE. .015) THEN 

SPEC1= 3.54272 - 522.1174 * AE + .1928462E+05 * AE**2 
ELSE IF (AE .LE. .025) THEN 

SPEC1= .3110063 - 25.2006 * AE + 520.0139 * AE **2 
ELSE 



- 80 - 



CS 

•K 

•K 

H 



cs 

•K 

•K 



tn 

fO 

fO 

a\ 



cs 



•K 

•K 

•K 

•K 

c 

H 

< 

Q 

»4 



o 

H 

2 

<! 

O 

O 

o 

Pd 

o 

Cz-( 

Q 

W 



w 


























W 








Pi 


































m 








CO 








M 


































• 








2 








2 


-/-N 




































m 




2 








O' 


Pd 




































• 




Q 








Pi 3 


2 


o 


























2 




2 






2 








^ Dd 


\ 


2 




CS 








00 








CsJ 






PJ 












O 








Pd 


2 






•K 








cr» 








CS 






2 




2 


2 






O 








2 


2 


•K 




•K 








un 








o> 






H 




M 


H 










2 




\ M 


v—' 
































2 




2 




o 


o 


> 




2 Pd 




m 




o 








CO 








O 


2 








> 




M 




M 


Pa^ 




• 


H < 


M 


<T» 




Pd 








• 








H 


W 




CS 






CO 


2 




CO 


2 




o 


2 


H 


• 




c 














O 




32 








1 




> 




c 






o 


WkJ 


< 


Vw/ 










• 








m 


O 


H 




• 






• 






2 


•K 


CO 


< 3 - 


C 


Pd 






•JC 






CV| 


1 






cs 


O 






O 




m 


O' + 








00 


o 


Q 2 












2 






















• 


w 






Pn 


2 


CS CO 


< O 


3 : 






2 


/<-s 


2 


•K 


cn 














• 






• 


m 




O 


H 


CN 


CS 


O M 


O 


H 


• 


M 


H 


< 




vD 








o 




to 




in 






• 






CO 


CS 




2 H 


2 


o^<r 


2 


2 


2 


2 


1— 1 










• 


• 




• 


•K 


<— ( 






2 


o o 




O 


Pa^ 




CS 


s.-/ 


s-x 


2 


< 


fO 






II 


• 


O' 










Ph 






O 




o 




H W 




+ 






H 






o 








H 


W 




O 




X 


o 


•K 




M 


•K 


• 


2 


< CO 


Pd 






\ 


2 


•K 


+ 


fO 








O 


• 




o 


•K 


< 


O 






H 






> 


2 1 


2 


M 






M 






o 






o 


• 










2 




X 






Q 


•K 




2 CO 


H 


o 


M 




2 


H 


2 


• 






o 




t-H 


2 


>»-✓ 


X 


> 


'wO' 


< 




2 


kJ 




•K 


CO 


<d 


H 


2 






2 


Q 






o 


kJ 


w 


o, 


M 


Cn 


<! 






2 




2 


2 


2 




Pd O 


:s 














II 








H 


o 


2 


1— 1 


2 


II 


M 


> 




O 


H 


2 


13 


W od 


















M 


II 


o 


O 


o 


> 




> 






s—/ 




M 2 


CO 


O 


<! 


CO o 


o 


II 


II 










o 




w 


CO 


2 






w 






w 




2 


< 








2 


2 


H 








II 




w 


o 


H 




S-/ 


S-X 


II 


CO 


II 


M 


CO 


II 


CO 


Q O 


ti 


II 


II 


Pi 3 2 


M 


2 


o 


li 


li 


H 


li 




2 


O 


o 






2 


2 


2 


Q 


2 


2 


2 


2 


H 


:s o 


Q C 


2 


o 


pd 


H 


H 


2 


2 


CO 


W 


2 


Q 


M 


H-l 


> 


M 


> 


> 


W 


> 


2 


Ci 3 -K 


O < 


:s 


2 H 


o 


H 


< 


2 


2 


H 


Q 




































O O 


o 








































■K 








O H 


o 








































•K 








II II 


II 










































cs 






H 2 


o 








































•K 


CS 






o< 


2 










































<y\ 
























o 
























CJ CJ O 








o o 


o 


O 













0D= DW + AE * 4 . 

DTEMP= (TSAT - TCI) / (TSAT - TOUT) 

0DTLM= (TOUT - TCI) / LOG(DTEMP) 

0TCS= (TOUT + TCI ) / 2 . 

CNU= .122181 - .148162E-02 * OTCS + .5516445E-05 * 0TCS**2 
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IF(INDEX1 .LE. 1) THEN 

LTOT= (9.125 + TIS) / 6. + TLT(KT) 

ELSE 

LTOT= TIS / 6. + TLT(KT) 

END IF 



TLT = TOTAL EFFECTIVE TUBE LENGTH ( FT . ) , ( TLT=TOTL ) 
LTOT= TOTAL TUBE LENGTH (FT.) 

LANE= STEAM LANE BREADTH (FT.) 
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END IF 

IF(C0NFIG .GT. 1) THEN 

WC0VER= .20A15 * DTB**2 * TCOVER * COVERD 
WHDR= WHDR/2 . + WCOVER 
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* (WCOVER + WUTUBE) + (LHW + LTOT/2 . ) * WTB) / WDRY 
ELSE 

DRYCG= (ARGIA + (LHW + LTOT/2.) * (WTB + WHDR)) / WDRY 
END IF 



ARG15= STMLD / 3690. * 61.5 
ARG16= WLIQ - ARG15 

WETCG= (DRYCG * WDRY + ARG16 * (LHW + LTOT/2.) + ARG15 * VHW / 

1 25.1328) / WET 

DRYCG=HEIGHT OF CENTER OF GRAVITY ABOVE CONDENSER BOTTOM FOR CONDENSER 
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TABLE 


4- 


Properties 


of Multiple 


-Helix internal 


Ridged Tubes 


Symbol 


Tube 

no. 


Nominal 
diameter, in 


Root 

diameter 

in 


Fin 

count, 

fins/in 


Outside 

area, 

ftVft 


ID 

(max)f//, 

in 


No. of 
starts 


A 


12 


3 


0.623 


(Stripped) 


0.163 


0.573 


5 


> 


30 


3 


0.627 


26.4 


0.640 


0.575 


6 


V 


22 


3 


0.625 


26 


0.640 


0.569 


6 


<3 


27 


3 


0.622 


26.5 


0.640 


0.572 


6 


A 


28 


3 


0.624 


26 


0.640 


0.573 


5 


► 


31 


3 


0.624 


26.6 


0.640 


0.576 


6 


T . 


29 


3 


0.625 


26.1 


0.640 


0.575 


5 


M 


38 


3 


0.633 


38.0 


0.330 


0.572 


8 


o . 


37 


3 


0.624 


38.5 


0.901 


0.574 


10 




40 


3 


0.627 


27.3 


0.689 


0.561 


10 . 


o 


41 


3 


0.628 


38.1 


0.852 


0.572 


12 * 


<3 


9 


3 


0.628 


38.5 


0.901 


0.575 


10 


□ 


21 


i 


0.740 


26 


0.640 


0.684 


6 


B 


19 


i 


0.745 


(Stripped) 


0.195 


0.692 


6 


/ 


44 


3 


0.627 


41.0 




0.574 


6 


\ 


43 


3 . 


0.628 


41.0 




0.573 


10 


o < 


42 


3 


0.626 


41.0 


0.901 


0.573 


10 


1 


46 


3 


0.627 


41.2 




0.573 


10 


( 


45 


3 


0.627 


41.2 




0.577 


10 


o 


13 


1 


0.883 


26 


0.841 


0.820 


6 


d 


32 


1 


0.877 


26 


0.871 


0.825 


6 


O 


25 


1 


0.880 


26 


0.841 


0.816 


6 


o 


24 


1 


0.878 


26 


0.841 


0.814 


6 


o 


23 


1 


0.880 


26 


0.841 


0.816 


6 


o 


26 


1 


0.864 


27.3 


0.841 


0.815 


6 



Symbol 


Internal ridging 


Internal aspect ratios 


ID of 

• envelope 
tube, in 


Height e, 
in 


Pitch p, 
in 


eldi 


^!p 


pMi 


A 


0.0125 


QMS 


0.0218 


0.0263 


0.829 


1.000 


> 


0.0162 


0.279 


0.0282 


0.0581 


0.485 


1.000 


V 


0.0165 


0.320 


0.0288 


0.0516 


0.558 


1.000 


<3 


0.017 


0.385 


0.0297 


0.0442 


0.673 


1.000 


A 


0.019 


0.475 


0.0332 


^ 0.0400 


0.829 


1.000 




0.0198 


0.287 


0.0344 


0.0690 


0.498 


1.000 


▼ 


0.0207 


0.469 


0.0360 


0.0441 


0.816 


1.000 


< 


0.018 


0.212 


0.0315 


0.0849 


0.371 


1.000 


o 


0.0200 


0.166 


0.0348 


0.1205 


0.289 


1.000 




0.0175 


0.170 


0.0312 


0.1029 


0.303 


1 .000 


o 


0.0215 


0.138 


0.0376 


0.1558 


0.241 


1. 000 


<1 


0.0204 


0.191 


0.0355 


0.1070 


0.332 


1.000 


□ 


0.017 


0.391 


0.0249 


0.0435 


0.572 


1.2348 


51 


0.017 


0.285 


0.0246 


0.0596 


0.412 


1.000 




( 0.021 


0.207 


0.0366 


0.101 


0.361 


— 




\ 0.021 


0.124 


0.0366 


0.169 


0.216 


— 


o 


< 0.024 » 


0.0949 


0.0419 


0.253 


0.166 


1.000 




) 0.021 


0.094 


0.0367 


0.223 


0.164 


— 




( 0.015 


0.094 


0.0260 


0.1596 


0.163 


— 


o 


0.0178 


0.333 


0.0217 


0.0535 


0.406 


1.5936 


o 


0.0193 


0.335 


0.0234 


0.0576 


0.406 


1.5936 


o 


0.0205 


0.340 


0.0253 


0.0603 


0.420 


1.5936 


o 


0.0205 


0.330 


0.0252 


0.0621 


0.405 


1,5936 


Q 


0.0205 


0.338 


0.0251 


0.0607 


0.414 


1.5936 


o 


0.021 


0.340 


0.0258 


0.0618 


0.417 


1.5936 
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TABLE 5 Friction Factor Characteristics 
of Multiple-Helix Internal Ridged Tubes 



Tube 



no . 


m 


r 


12 


0.762 


0 


30 


0.64 


- 0.00039 


22 


0.697 


- 0.00017 


27 


0.72 


- 0.00026 


28 


0.72 


- 0.00028 


31 


0.61 


- 0.00095 


29 


0.68 


- 0.00064 


38 


0.61 


- 0.00109 


37 


0.54 


- 0.00392 


40 


0.627 


- 0.00080 


41 


0.57 


- 0.00259 


9 


0.59 


' - 0,00197 


21 


0.70 


- 0.00014 


19 


0.626 


0 


44 


0.53 


- 0.00180 


43 


0.55 


+ 0.00017 


42 


0.58 


+ 0.00750 


46 


0.54 


+ 0.00275 


45 


0,52 


- 0.00159 


13 


0.63 


+ 0.00024 


32 


0.70 


+ 0.00090 


25 


0.645 


+ 0.00032 


24 


0.64 


- 0.00035 


23 


0.637 


+ 0.00028 


26 


0.64 


+ 0.00018 




FIGURE 17 Friction Factor Curves for 
Multiple -Helix Internal Ridged Tubes 
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FIGURE 18 a 

Operands for friction factor equation versus geometric aspect 
ratio e/l for multiple -helix internal ridging with p/d^> 0.36 




FIGURE 18b 

Operand for friction factor equation versus geometric aspect 
ratio e/l for multiple-helix internal ridging with p/d^< O.36 
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WATERBOX AND TUBE END LOSSES 
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